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ABSTRACT 


The dissipation method is uSed to obtain estimates for 
the friction velocity, Uys, as well as values for the neutral 
drag coefficient, Cp, for data collected from a coastam 
tower off San Diego, California. Cpy is found to be indepen- 
dent of the ten-meter height windspeed, Uj9, for velocities 
between 4-9 # m/sec. Its value ius estimated 1c) be 
(0.94+0.4)102 which compares well with values by Smith 
(1980) and Barge andePond( 1931) Pettite  tnendomin eG , 
with fetch and sea state are also observed. Drag coeffcient 
estimates are found to be higher for short fetch than for 
long fetch conditions. Cp, is also seen to increase sharply 
just before frontal passages and during sea breeze condi- 
tions when the waves are actively growing. With the wind- 
speed and wave fi ec eeete rine [el tele ye Pee CT yay hs et Ol ica 


decrease with time to a smaller and more constant value. 
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I. INTRODUCTION 


The Marine Atmospheric Boundary Layer (MABL) is the 
region adjacent to the air-sea interface where the vertical 
variation of the mean stresses is large compared to the 
horizontal gradients (Holton, 1979). In this layer, turbu- 
lent shear stresses due to both molecular properties 
(viscosity and diffusion) and turbulent eddies coexist. 
These eddies are very effective mixing agents which serve to 
transfer momentum, heat and other quantities to and from the 
earth's surface, at a rate much faster than that for molec- 
ular diffusion. Frictional forces due to molecular viscosity 
are usually neglected by means of scale analysis (Holton, 
1979), and the emphasis has been placed on parameterizing 
the fluxes resulting from the turbulent eddies. 

The mean vertical fluxes of momentum, sensible heat and 


latent heat are (Fleagle and Businger, 1980): 





T = -pgu Ww! (la) 
H = C40 ww’ Gey) 
Jd) = Lew’ q' Gle) 


where 7, H and E are the momentum flux (wind stress), 
sensible heat flux and latent heat flux, respectively. u' 
and w' are the horizontal and vertical turbulent wind compo- 
nents, while @' and q' represent the turbulent fluctuations 
of potential temperature and specific humidity. ie ts) ene 
Paeeme meal Of vaporization, p, is the air density and Cy, is 
the specific heat of air at constant pressure. 

In order to estimate these fluxes using routine surface 
layer observations, bulk transfer formulae were developed 


(see Fleagle and Businger, 1980): 


T is palp (Uz = ss (2a) 


dial 


H = meee, 1 (hy = Be) (2b) 
= PalsCR(Uz - Ug) (dg - dz) (2c) 


where U,, T, and q, represent the windspeed, temperature and 
humidity specified, by convention, for a reference height, 
z- U,, T, and q, are the same quantities measured at the 
surface. Cp, Cy and Cp are dimensionless quantities repre- 
Senting bulk transfer coefficients for momentum, sensible 
heat and latent heat, and are for the height, z. These are 
known as the drag coefficient, the Stanton Number, and the 
Dalton Number, respectively. 

The turbulent exchange processes greatly influence the 
general circulation of the atmosphere and thus have tremen- 
dous impact on countless social and environmental concerns. 
Central to their understanding is a determination of the 
nature of the bulk transfer coefficientsyparticulgmy the 
drag coefficient, Cp. This quantity encompasses the mecha- 
nism responsible for the wind stress at the surface. Since 
most of the parameters used to describe the bulk formulae 
are routinely measured, ascertaining a@uuetiable values ouse, 
based on these quantities is essential for accurately 
defining the surface stress, 7. Further, Cp is an important 
factor driving many of the current ocean wave, acoustic and 
boundary layer models (Geernaert, 1985a). 

Considerable effort has gone into examining the over- 
water drag coefficient. It has been shown to be windspeed 
and stability dependent, and may also vary with factors such 
as fetch and depth (Geernaert, 1983). Studies also indicate 
that the overwater drag coefficient tends to have a magni- 
tude on the order of 107°. Its value exhibits considerable 
Scatter under varying environmental conditions. Higher 
values are usually associated with an unstable atmosphere, 
short fetch and/or shallow seas, implying a wave dependence 


(Geernaert, 1983). The more recent use of remote sensing to 


Wye 


Study winds over the ocean requires an understanding of this 
wave dependence since it assumes that the wave structure is 
related to the wind stress. Clearly, the processes 
affecting the drag coefficient merit considerable further 
study. This thesis will investigate the dependence of the 
drag coefficient on windspeed, fetch and wave state. 

An important scaling parameter, based on the wind 
stress, is the friction velocity , Us, which may be defined 
as: 


T = 1 ples ( 3 ) 


The quantity Us is related to most air-sea exchange 
processes such as heat and moisture fluxes (Businger, 1973), 
the surface drift (Hicks, 1972) and ocean mixed layer 
dymamacs (Kraus, 1972). 

Data considered in this investigation were from the 
Tower Ocean Wave and Radar Dependence (TOWARD) experiment. 
Located offshore of Mission Beach at San Diego, California, 
the Naval Ocean Systems Center (NOSC) tower provided a 
Stable platform with a capability for mounting a variety of 
sensors to measure the variables needed to parameterize the 
bulk formulae. Naval Postgraduate School (NPS) personnel 
made atmospheric surface layer measurements of both mean and 
fluctuating velocity, temperature and humidity. The NPS 
measurements are listed in Table I. Wave information based 
on radar Signal returns was collected through investigations 
from the Naval Research Laboratory, the Jet Propulsion 
Laboratory and the University of Kansas. The Scripps 
Institute of Oceanography measured tower frequency wave 
Spectra with an array of sub-surface pressure transducers. 

The measurement strategy for the TOWARD Experiment was 
originally composed of two intensive data collection 
periods: September/October 1984 and January/ February 1985. 


These periods were chosen because of their differing average 


LS) 


windspeeds and resultant sea states. With its relatively low 
average windspeeds, the fall measurement period had calm to 
moderate sea state conditions. The winter period, with its 
higher average windspeeds, was designed to focus on more 
active sea states. Equipment problems, however, forced the 
postponement of Phase I to October/November 1984 and Phase 
il to Maren) Aprivaeicco- 

Figs. 1.1 and 1.2 provide two views of the NOSC tower. 
Located in 17-18 meters water depth, the tower provided a 
Stable platform for data collection. A structural extension 
of the tower from the south side allowed wave rider sensors 
to be placed away from the main tower structure. 

Measurements of wind speed and direction, wave height, 
tide level and near surface currents were recorded routinely 
during intenSive measurement periods. Air and sea-surface 
temperatures were also recorded, along with wind stress data 
that was obtained with a sonic anemometer and a hot film 


S@EMNnSOTr. 
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TABLE I 


NPS Measurements 


Measurement Sensor 
Sea-Surface Temperature Floating platinum 
thermometer 


Mean Surface Layer: 


Wind (speed and direction) Cup anemometer, 
vane 
Temperature Aspirated platinum 
a: thermometer 
Humidity Dew cell | 
SS cooled mirror) 
Radiation (solar and IR) ppleys_(up_ and 
down), Total 


Turbulent Surface Layer: 


Wind siteress: direct Sonic anemometer 
Wind stress: indirect Hot films 
| En oLsSipation rates Hote films 
Sensible heat flux Platinum thermometer 
and sonic anemometer 
Latent heat flux Lyman-alpha and 
Senremenemonecer. 
Spectra: windspeed sonic anemometer 
temperature Platinum thermometer '| 
humidity Lyman-alpha | 
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Figure 1.l The NOSC Tower. 
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Figure 1.2 The NOSC Tower. 
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II. REVIEW OF EXISTING SURFACE LAYER THEORY 


The surface layer comprises approximately the lowest ten 
per cent of the MABL and is defined as the layer of air 
nearest the ground where atmospheric parameters may be 
scaled with surface fluxes. Within this layer, an approxima- 
tion adopted is that the fluxes are independent of height 
(Bustnmgeme 19/3. | 


In Eqn. (2a), the surface speed, U is a parameter 


Oo b 
necessary for estimation of the drag coefficient. However, 
because U, is very small compared to Ujg, U, is assumed to 
be zero in the bulk definitions. With this assumption, Eqns. 


(2a) and (3) may be combined into the form: 
[p10 = (Ux/Uy9)? (4) 


Note that the drag coefficient in Eqn. (4) is defined for a 
hel ghee Ol meters anidmloeden@eed 42S Chip) 
The similarity theory which described the mean surface 


layer wind profile is (Monin and Obukhov, 1954): 
QU/dz = (Ux/kz)o_(2/L) ©) 


where QU /dz represents the vertical gradient of the wind 
speed and k is the von Karman constant, which ranges in 
value “from .35 to “43 (Panotsky and Dutton, 1934). Thats 
study assumes a value of 0.4 for this constant. ¢,(z/L) is 
a dimensionless stability function. One parameterization of 
EXiceerumection is that Spresented Voysebusineeunemrect a la 
CLOW ae 


dm = (l+=a(z/D) “Gro (Be) 
Om = 1) 282 (z/L>0) (6b) 
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where q@ and B are diabatic profile constants, with approxi- 
mate values of 1541 and 520.5, respectively (Geernaert, 
1983). 2 is the meaSurement height and L is the 
Monin-Obukhov length, defined as the height at which the 
buoyant energy production equals the shear production 
Pous@mger, 1973). z/L is defined as: 


z/L = -gkew'T,"/T,' (Ux?) (7) 
Sveevirvvals temperature Fluctuation, T,, 1s related to 
the air temperature fluctuation, T', and specific humidity 


feeetuation, q , by: 





a ee (iteolg) + .61Tq’ (8) 
Thus ; 

wet. eel (nolan + .61Tw q_ Co) 
where g is the gravitational acceleration. The first term 


on the right hand side of Eqn. (9) represents the tempera- 
ture flux and the second term describes the humidity flux. 


These terms are further estimated by bulk equations where: 





wT! = CyUi9(T, oy ey) (10a) 
t ' - CrUig(d, = Ste) (ib: ) 


= 
Ne. 


Both Cr and Cy are on the order of 1079, but are less well 
defined than the drag coefficient (Geernaert, 1983). 


Umdier neutral or adiabatic conditions, 
bn(2/L) = ail eri) 
and Eqn. (5) reduces to: 
QU/dz = Tees (2) 
imrerratinoebaon...{ 12): 
u{z) 


fey AU = (Ue/ky fe dele a3) 


yields the well known logarithmic wind profile: 


BES) 


U = (Ux/k)1n(2z/z,) (14) 


where the roughness length, z is a measure of the effec- 


0? 
tive surface roughness. 

For mnon-neutral (diabatic) conditions, the stability 
parameter must be included. Integrating Eqn. (5) yields the 
Stability dependent logarithmic wind profile (Paulson, 


S707) = 
U = (U,./k)[in(z/z,) as W ] (139) 


where yw represents the stability parameter, 


For unstable stratifications (z/L<07) 


Ww = 2in{(1+x)/2] + In[(1l+x2)/2] - 2(tan)7tx + (7/2) (16) 
where: 


x + ee alam (16a) 
For stable conditions (2/L=>0): 


y IGA) (17) 


The logarithmic wind profile iS iperbant primaries tom 
two reasons. First, it can be used to estimate the ten 
meter height wind speed, Uj;g- Beginning with Eqn. (14), the 


Wind speed at a particular height, z, may be described as: 


o = (Ux/k)[1n(z/z,) - 2] (18) 


At ten meters elevation: 


TG =  (U)/kP (10/22 ee ee Gk) 
Combining Eqns. (18) and (19) yields the following relation- 
Sip 

Uio = U, + (Ux/k) [nC 10) Zoe (20) 


THis particular equation iS 9amportantemmecauce=tL se cemmot 


always possible to place sensors on ships or platforms at 


20 


the required ten meter height. Hswea Cesult jee l)@e Ccanwebe 
estimated regardless of the height of the instrument. 

Second, the log wind profile is important in determining 
the neutral drag coefficient, Con» Which may be defined as 
the wind stress exerted on the surface for neutral stratifi- 
cations. Estimation of the neutral drag coefficient is 
important because it eliminates variations in Cp parameteri- 
@eeLons due ~to differing stability conditions. Many 
modellers prefer to be provided with Cpy instead of Cp 
because the in-situ bulk weather data may be used to esti- 
mate a stability correction. Eqn. (24) may then be applied 
Boma given Chy andey. 


A. GENERAL DRAG COEFFICIENT REPRESENTATIONS 
Foams. (4) and (15) yield: 
Cy = (Ux /U,)*(k/[1n(2/25) - p_(2/L)]]? (21) 


Ender Ggeeutral conditions, where y=0: 


ane ei, (zy 2a (22) 
ioe oa eonoec,, 42re thewdrag coefficient and roughness 
length under neutral conditions. Rearranging Eqn. (22) 
ylelds: 

Mee o) = iIn(z/z,,) (23) 


Combining Eqns. (21) and (23) produces the following general 
result (Geernaert, 1983): 


o Se ge lec Toyl(Zee/Z,) 1° (24) 


The drag coefficient is known to vary under differing 
environmental conditions. Stmdieés Hayvemandicated that Cy 
increases more rapidly with increasing winds over water than 
over land (Hsu, 1974). Geerneart (1983) discussed the 


apeemettetoumoceldqn. (24) to both land and water cases. Over 


a= 


land, the roughness length is assumed to be largely indepen- 


dent of the surface stress.” Thus, 2) -meljgemd Eqnime( 24 ae 
reduced to: 
Cy =" (CpyP/* - ps)? (25% 


Over water, both the drag coefficient and, hence, the 
roughness length show considerable windspeed and stress 
dependence (Kitaigorodskii, 1973). Chasnock (19 35emeewac 
among the first to recognize that Cp increases considerably 
with windspeed over the sea. He postulated that the rough- 
ness length, and hence, the neutral drag coefficient was 
proportional to the wind stress through the following 


scaling argument: 
z = q'Uy*/ 26 
O ms 04 we 24 . ( ) 


where qa is the Charnock constant with a typical value of 
.015-.020 (Geernaert, 1983). Charnock's work served as a 
basis for subsequent parameterizations of the drag coeffi- 
cient and the roughness length in terms of sea _ state 
elements. 

Kitaigorodskii (1973) further described z, based on wave 
Spectral densities. Since that time, 2, and Cp have been 
shown to be dependent on a variety of other factors that are 
linked to the character of surface waves such as fetch, 


depth and wave age. Cp, therefore, can be evaluated as: 
Cy = Cpy(Uj9, 2/L, waves) C27) 


As a result, a different relationship is often used to esti- 
mate the drag coefficient over the ocean. - Applying 


Charnock's relation to Eqn. (24) yields (Geernaert, 1983): 


Cp = [Cpyt/* - w/k + (1/k)In(Cpy/Cp)]7* = (28) 
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Be DETERMINING THE WIND STRESS BY THE DISSIPATION METHOD 


Since estimation of the drag coefficient requires values 
for the friction velocity, it is important to examine how 
the momentum flux is measured. Two primary methods for 
obtaining - are the eddy correlation technique and the 
dissipation method. The former method involves measurement 


t 


of the turbulent wind components u’, v' and w', and ideally 
requires a steady platform. Because of the low frequency 
motions inherent in sea-going vessels on which most marine 
Studies are conducted, the dissipation method is often 
substituted for oceanic wind stress experiments. 

The dissipation method is an indirect approach in which 
the rate of dissipation of turbulent kinetic energy 1s used 
to determine U». This method requires application of the 


Bau ent kinetic equation: 
@e/Ot = Ux*QU/dz+(qw'T,')/T,-0/dz[w'e+(l/p,)w'p'J-e = (29) 


VheGewge/otems the local time rate of change of turbulent 


kinetic energy (TKE), Ux QU/ dz is the wind shear produc- 





Exon, Cawest )/ 7. represents the buoyancy Romeese 
ogee eG oew Pp | ise the redistribution of the flux of 
TKE and work done by pressure fluctuations and ¢€ is epsilon, 
the rate of dissipation of turbulent kinetic energy. 

Steady state conditions are assumed. The divergence term 
is small relative to the remaining terms in the equation and 
Pere renorcmmece lected (Large and Pond, 1981). Eqn. (29) 


ehemesimpliires to: 


Ux-dU/dz + [(gwTy')/T,] - € = 0 (30) 


Saeoreecompbinime cams. (5), (7) and (30), and rear- 


ranging, an expression for the friction velocity is 
obtained: 
= 7 WS) (a) 
Ux = [ekz/{¢6,(2/L) - (2/L)]J] 


Zo 


For unstable conditions, the denominator is often replaced 


by: 

@-(2/L) = dyl2z/L) - (2/L) (32) 
Eqn. (31) may then be written as: 

Hie - [ekz/$-(2/L)]*/° (33) 


with stability corrections for dissipation scaling (Wyngaard 
and Gote: nto se 


be = (1+ 0.5}z/L}3/2)3/2 = z/L<o (34a) 
$< = (1 4225) 2yi-' =)! = ergo (34b) 
The friction velocity, Us, can therefore be estimated in 


terms of the dissipation rate, €. Kolmogoroff theory is then 
employed to determine a value for e€ based on high frequency 
wind fluctuations. All spectra can be divided into three 
portions: the energy producing subrange, the inertial 
subrange and the dissipation subrange (Panofsky and Dutton, 
1984). Large eddies from the energy containing subrange 
cascade kinetic energy and momentum down to smaller eddies, 
eventually reaching the dissipation region where molecular 
viscosity converts the kinetic energy to heat. Between these 
two regions is the inertial subrange. In this region, energy 
is neither produced nor dissipated. Kolmogoroff theory 
States that in the inertial subrange of isotropic turbu- 
lence, the high frequency fluctuations of the wind velocity, 
u' can be expressed by a one-dimensional power spectrum 
(Sechacher, et. ai tle. 


Su (k., ) = awet! Ik, 9/3 Coy 


where Su is the energy density spectrum of the horizontal 
Windspeed, is a coefficient of magnitude 0.50 (Champagne, 


et. al., 1977), and ks is thevwavenunmbpen 
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Applying Taylor's frozen turbulence hypothesis: 

Ke =  Qqrf/U (36) 
equation (35) becomes: 

Su (k,,) = axe?! 3 (2 £/B) 77/3 ian) 


The Kolmogoroff Spectrum can now be integrated between 
two wavenumbers, designated Kxe1 and Ken to give the energy, 


o- Aes mtewval(Schacher, et. al., 1981): 


Kora 
ae Bie Wear dic. (38) 
= imoqeee Gee, 6/2 = fen) (38a) 


comemine Bams. (36) and (38) yields: 
oA gE = 0.230(e0)2/3(£,-2/3 - £,°2/3) (39) 


where f£4=(ky1U)/2 and f= (kepU)/2. Red@mrancine Equa (39) by 
Heme Eg. (35) produces: 


Ux = 2.65oyqclz/ lle. (2/Ly)Jt/?(£ °°! - £,°°/ 2) 1/4 (40) 


Peaeetcemnew et De experiment, £7; = 5 Hz and £, = 50 Hz, and 
z=22 meters. Substituting these values reduces Eqn. (40) to: 


ie = 7.8loya slUg_(2/L)]~*/° (41) 


Because o_(2/L) also depends on Us, Eqn. (41) must be 


solved iteratively. 


Cee SLU HOT FILM CALIBRATION 


Calibration of the hot film sensor was accomplished by 
recording the changes in voltage produced by the winds 
passing over the sensor. The relationship may be expressed 


ba 
ve = ge + BUrell/4 (42) 


JES) 


where V is the voltage, V. is the plotted calibration inter- 


cept and Urel is the eclaeeen windspeed. B represents a 
calibration factor that must be continuously determined for 
each data collection period. 

B may be esStimated either through the dynamic or the 
regression methods. The dynamic method begins by differen- 


tiating Eqn. (42): 
B = 4vUrel!/2(dv/dv) (43) 


Over each data collection period, dV/dU is estimated by 
ovy/og> where on and Oy are the standard deviations of U and 
V, respectively. Averaging V and Ureleand solvine tors 


mae lds: 
B = 4vUrel+!* (gy/o7) (44) 


Rearranging Eqn. (44) and substituting into Eqn. (40) 


produces the following estimate for VU»: 


Uy. = 
2.81[(4V/B)Urelt/? Joy, glz/Up_(2/L)] 1/2 (£7277 ~~ ft peaaeD 


The regression method is based on values of V and Urel 
over several time periods with varying average relative 
windspeeds. V is then usually plotted as a function of Urel. 
According to Eqn. (42), the points showla tall alene@a lame 
with a slope, B, where Voy Sepresemtrs the y-intercept. 

Measurement-relative scatter in both Urel and V can 
cause errors in the slope. Therefore, it 1S Sometimes advan- 
tageous to choose V as the independent variable if it shows 
less scatter. In cases where it is difficult to determine 
which variable is the least reliable, two regressions are 
calculated, with Urel and V alternated as independent vari- 
ables. B is then determined from the average of the two 


slopes. 
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D. WIND-WAVE COUPLING DYNAMICS 


The total wave field energy is divided into swell and 


wind wave components: 


a Ge ig (46) 
where oT represents the total energy in the wave field 
and of. and Ong are the energies contained by the swell and 


wind waves, respectively. o2 


2 


Po cConmolmpe predicted bacedacn 
local scaling. o°,,, has received considerable study and a 
great deal of theory exists concerning this parameter. The 


approach currently is to describe a2 


ee Ens) OL ideal 
equilibrium spectrum. This theory assumes that wind forcing 
controls the dynamics of every wave component of the wind 
wave spectrum. Each wave is characterized by an angular 


frequency, n’', where: 
n! =a £ (47) 


Gravity is the primary restoring force for theSe waves. 
The local equilibrium energy density spectrum, S(n'), is 


defined as: 


ae = a Sims sen” (48) 


and S(n') may be parameterized in terms of Phillips scaling 
arguments (Phillips, 1980): 


S(n') = Bg(n')~? n'>n,' (49) 


Beta is the Phillips coefficient which can vary from .006 to 
fee Gee materia, 91703). mn, represents both the wind-wave 
peak frequency and the boundary between the swell and equi- 
librium ranges. Waves of frequency ES | represent swell. 
Wind waves are defined to be those with frequencies of 


ie, . 9, te usually approximated by: 


' 


Za 


for steady state conditions (Kraus, 1972). 
Beta is often described in terms of dimensionless fetch, 


ra 
5 =  gx/ (Uy) (51) 


where x is the upwind fetch. Large values of B are associ- 
ated with small values of fetch. Beta has also been parame- 
terized over an applicable range of .25 to 1.0 hz by 
(Geernaert, et.al., 1985a): 


B = 0.005 + .002n’ + 1.5(C,/Uy.) 72 (52) 


where C, represents the phase speed of the largest wind 
wave. (C,/Ux) is defined as the wave saturation parameter, 
and this parameter behaves in the same manner as the wave 
age C/U. (C,/Us) is a dimensionless indicator of the degree 
to which the wind waves have reached steady state and is 
frequently related to the non-dimensional fetch as (Wu, 
JSkeisy 


(Gh/ Uy = y0n05 Cex nC sates: (53) 


Hsu (1974) related the wave saturation parameter to the 


roughness length, 2z by defining the Charnock coefficient 


O°? 
in terms of the slope of the wave of frequency ne = His 
formulation became: 
Zo = (H/2qr)(C,/Ux) 77 (54) 
where H is the significant wave height. Z~ can be seen Jise 


O 
vary according to the local wind wave spectral peak. 


Numerous wave models have incorporated these relation- 
ships in an attempt to relate variations in the surface wave 
energy spectrum to variations in the magnitude of the drag 
coefficient. Davidson (1974) formulated an empirical model 


in which the drag coefficient is represented by: 


Cp = ([k/[1n(z/2,) - p + b{(Co/Ux) - 26.3]]]? (55) 


Zo 


where b is an empirical constant found to be 0.13. 
Others have employed indirect methods to calculate Cy 

uSing roughness length parameterizations. Byrne (1982) 

derived a roughness length in the form of an integrated 


Slope spectrum: 
o t t t 
Zo = Caf S(n')(n')@adn (56) 


where Goes emeaeeusteb Wes coefficient . 
Kitaigorodskii (1973) derived a roughness length model 


in the following form: 


i = 2/30(f™ S(n' Je" 2k8/n'Usdn' )1/2 (57) 


One purpose of the TOWARD Experiment was to evaluate 
Scatterometer return Signals in order to examine surface 
roughness as a predictor of windspeed over the ocean. These 
Signals can be correlated with windspeed since they increase 
with increasing surface roughness. The surface roughness is 
dependent on tyewhich from»Eqn. (2a) is a function of Cy and 
the ten meter height windspeed, Ui eee Tesicat Leronetw et, 
however, evaluates Ujo based on signal reflection from the 
Surface (Geernaert, 1983), while the drag coefficient has 
been treated as a constant or as a simple function of Uig- 
Thus, the dependence of the drag coefficient on the surface 
roughness is not fully understood. Because wave growth and 
deecy waiGiect thewmagnitude ofeCp (Geernaert, 1985a), Sscat- 
terometer predictions of windspeed based on the bulk 
formulae could be improved by eh nacne roughness as an 
additional variable. 

The TOWARD Experiment was designed, in part, to look for 
Geecemdein Cy 4s agfunction of time over the course of the 
sea breeze. The sea breeze is characterized by growing winds 
and sea state along with corresponding increases in surface 
roughness. By observing the behavior of this coefficient 
With time, the interaction between wind and wave state may 


be better understood. 
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LII. SYNOPTIC METEOROLOGY DURING THE TOWARD EXPERIMENE 


The TOWARD experiment was conducted at the Naval Ocean 
Systems Center (NOSC) tower which is located two kilometers 
offshore of Mission Beach, San Diego. Minimal surface layer 
climatological data is available for the site. Blanc (1981) 
summarized the environmental conditions for San Nicholas 
Island (SNI), California, which is loeated rouehty  '50 ian 
northwest of the NOSC tower. While wind and wave conditions 
are more vigorous at SNI than those at the tower, the air 
and sea-surface temperatures are consistent between the two 
Sites. Blanc determined the most frequent true wind direc: 
tion to be from the northwest throughout the year with mean 
windspeeds of 7.7 m/s in March and 6.2 m/s .for both October 
and November. During the months of March, October and 
November, when the TOWARD experiment was conducted, the mean 
climatological daily temperatures (9C) are 15.8 for March, 
20.4 for October and 18.4 for November, with average daily 
minimum temperatures of 12.4, 16.0 and 14.9, respectively. 
Mean climatological sea-surface temperatures for the same 
months are 14.4 for March, with values worms .2 and. l6. 6G jtem 
October and November. 

The mean conditions for San Nicholas Island differ some- 
what from those reported onshore at Lindbergh Field in San 
Diego. The mean daily maximum temperatures for March, 
October and November are 18.9, 23.3 and 21.2, respectively, 
while the average daily minimum temperatures are 10.1, 14.7 
euovol INO/es) - In March and October, the prevailing winds are 
generally west-northwesterly but shift to a northeasterly 
direction in November. Mean windspeeds are also smaller 
with values of 3.3, 2.9 and 275) Ssecmreported For rene 


October and November, respectively. 
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The TOWARD Experiment was designed to cover a large 
Pomee SOL Wwindee conditions, particularly for high winds. 
Measurement periods were therefore chosen to be from 
mid-September through mid-November 1984 (Phase IL) and March 
1985 (Phase IIL). 

The almost daily presence of the local sea breeze is a 
factor in the environmental conditions at the NOSC tower. 
Sea breeze circulations develop when air over land is heated 
much more than the air over the ocean. Since the land 
usually becomes warmer than the sea during the day, density 
differences over the two surfaces bring cool air from the 
sea toward the warmer land at low levels. A return flow of 
warm air from the land to the ocean occurs at upper levels. 
At night, radiative cooling takes place over the land while 
the sea surface maintains a more uniform temperature due to 
the large heat capacity of the ocean. The result is 
nocturnal, low level flow from the land toward the sea 
(ess, 1979). 

Because the largest horizontal density gradient occurs 
at the coastline, the sea breeze dominates in the coastal 
zone, while the geostrophic wind (ve) becomes a more impor- 
tant factor for describing surface layer energetics farther 
offshore. It follows that waves generated by the sea breeZe 
wouwtlantend so be of relatively short fetch. Fig. 3.11 from 
Seal ceober serovides an example of synoptic conditions for 
Southern California during the TOWARD Experiment. The large 
scale flow (V,) can be seen to originate from the northwest. 
Smepmewopservations, 200-400 km offshore, also indicate a 
northwesterly flow. These wind observations from offshore 
Seees Gchoweslattle Saiurnal variation. At 1000 local time, 
however, reports from Lindbergh Field, San Diego, show winds 
from a south-southwesterly direction, indicating smaller 
scale, local-forcing nearer the coast. A significant diurnal 


shift in wind direction can also be seen in theSe reports. 


on 


Thus, the sea breeze appears to exert an influence less than 
300 to 400 km from the shore. 
Since the sea breeze is of a different scale stnan te 


Synoptic, geostrophic wind, computation of V, provides an 


= 
indication of the relative magnitudes of these two types of 
forcing near the NOSC tower. Calculations of the geos- 
trophic wind from synoptic charts were made for 1000 local 


(1800 GMT) and 1600 local (2400 GMT). 


Ve =  (1/£4) (dp/d®) (58) 


kh 


where Ve represents the geostrophic wind, f is the coriolis 
parameter, p is the atmoSpheric pressure and n is a unit 
vector normal to the isobars. These times were chosen in 
order to compare the differences in sea breeze magnitudes 
between morning and afternoon. Typically, the strongest sea 
breeze occurs in mid-afternoon when the land surface has 
heated to its daily maximum. Since #f£imal sunfacemechagaiae 
were not available for 2100 GMT { 130G@Miocal).. analvyoesc sion 
2400 GMT were uSed instead. These values were compared with 
windspeed observations collected at the tower and with data 
from ship reports in the Southern California area. 

Data collected at the tower indicated a general increase 
in windspeed from 1000 to 1600 local (1800 GMT to 2400 GMT). 
The geostrophic wind velocities based on Eqn. (58) and the 
Synoptic pressure maps were smaller and showed much less 
diurnal variation than those of the tower CORED Ve Ship obser- 
vations from 300-400 km offshore indicated consistently 
higher windspeeds for both time periods, usually ranging 
from 5-15 kt. Although the computed Vo 


surface winds (V,) showed considerable differences in magni- 


and ship reported 
tude, the wind directions from offshore sites during the 


experiment showed no statistically significant diurnal 


effects. Tables [II and LIL showeemnese Jresutes. 


a2 


TABLE II 
Computed Ve & Ship Reported Sfc Winds (Phase I) 
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TABLE III 

ak, 
Computed Ve & Ship Reported Sfc Winds (Phase ITI) 
(V+; Ve and V; in m/s) 
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A. THE TOWARD SYNOPTIC AND SURFACE LAYER CONDITIONS 


This section will be a presentation of the surface layer 
and larger scale meteorological conditions for the TOWARD 
experiment. The data set is divided into two sections, Phase 
I and Phase II. Phase I covers the period from 2 October to 
8 November 1984, while Phase II covers the period of 12 to 
Mace Wlelveela, IRS), Each phase is discussed separately. 


Individual surface analyses for 1800 GMT will be presented 
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Peene  tecemmLOt “days On Which data were collected. 
Additional figures representing wind speed and direction, as 
well as air and sea-surface temperatures can be found in 
Appendix A. -The latter figures present time series of 
discussed quantities and provide a reference for’ the 
following discussion. Unless otherwise noted, the times are 
local (PDT for Phase I and PST for Phase II). 


1. Phase I 


During the first week of the TOWARD experiment, 
winds were from the northwest and windspeeds were light to 
moderate. A maximum value of roughly 7 m/s (Fig. A.1) was 
recorded on October 6th. Strong local diurnal variations due 
to the sea breeze also prevailed in the first week, which is 
meee rcularily veyident in the data from 3 and 6 October. On 
those dates, windspeeds increased Significantly from 0900 
PDT, remained steady for several hours, and then decreased 
markedly throughout the sampling period which ended at 2215. 

Data collection resumed at 1235 PDT on October 6th. 
Recorded windspeeds remained fairly steady from this time to 
1500, when values began decreasing. Data were only collected 
for one hour on 7 October. Windspeeds increased steadily 
from about 4 to 6 m/s over this period. 

S¥Mapelemenarts forethe 2-8 October period do not 
indicate any significant weather activity. A low pressure 
center was present southeast of the San Diego area early in 
the week, before moving farther southeast. Weather condi- 
tions were subsequently dominated by a strong offshore ridge 
which persisted throughout the week. 

Data were collected only on the llth and 12th of 
October during the second week of TOWARD. Wind directions 
were from the northwest on the llth and windspeeds were low. 
On the 12th, winds were also out of the northwest, but had 


much higher velocities. A maximum of about 9 m/s_ was 
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Figure 3.1 Surface analyses Taco GmMr 2 (Cay 
and 3 (b) October 1984. 
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BRegure 3.2 Surface analyses 1800 GMT 6 (a) 
and 7 (b) October 1984. 
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recorded on that date at 1550. Typical sea breeze patterns 
were also apparent in the data partremtariyeeon sene 2th 
Weak easterly winds were present in the early morning. The 
winds shifted to the northwest near noon and increased 
rapidly to a maximum. These values began to drop off at 
approximately 1600 PDT. 

The synoptic charts from this period show a frontal 
passage through the San Diego area on the morning of 12 
October. This event cannot be correlated with the strong 
winds of the 12th, because the front was already several 
hundred kilometers to the southeast when the winds were 
recorded. Synoptic charts do not indicate the presence of a 
Strong, large-scale pressure gradient associated with the 
trough. Thus the most likely cause of these high winds was 
Strong diurnal heating and subsidence following the passage 
of the frome. 

Data collection for the third week of the experiment 
began at Q900 PDT on October 15th. A definite wind shift can 
be seen on the 15th accompanied by sharp fluctuations in 
windspeed. Synoptic charts for that date indicated the pres- 
ence of a stationary front almost directly over the San 
Diego area at 1000 PDT. This front passed through the area 
at a rate of only 2 m/s. Its slow rate of movement caused it 
to influence weather conditions in the San Diego area 
through the 16th. 

Wind shifts plus strong increases and decreases in 
windspeed were also apparent on the 16th and 1/7th. Synoptic 
charts for these dates show a front moving rapidly through 
the tower area at a rate of over 10 m/s. This front was well 
to the southeast within six hours: 

October 18th provided a more typical weather pattern 
with northwesterly winds and a strong increase in windspeed 
between 1000 and 1700 local due to the diurnal heating 


pattern. Synoptic charts show a strong high pressure ridge 


38 


/ 


SMe] CS Lt. 
Tae Hes 


Py: 


Fee RS 


= = 





Figure 3.3 urface analyses cere GMT 11 (a) 
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Figure 3.4 Surface analyses 1800 GMT 15 (a) 
and 16 (b) October 1984. 
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Surface analyses 1800 GMT 17 (a) 
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and 18 (b) October 1984. 


Figure 3.5 


building to the southwest, with a weaker ridge to the 
northeast. 

Winds were westerly on the morning of the 20th with 
a shift to an easterly flow near midnight. Velocities were 
moderate and fairly constant, with values decreasing in late 
afternoon. 

The first observations for the week of 22-29 October 
began at L1OO on October 23rd. A strong increase in wind- 
Speed can be seen, along with a northwesterly wind direc- 
tion. On the afternoon of the 25th, windspeeds again showed 
an increase. The wind direction was southerly in the late 
morning and shifted to the northwest by early afternoon. 
Wind speeds recorded on the 2/7th were moderate with little 
diurnal chamee. This pattepm@ecentinucdsenmencez ema, althoued 
the winds again originated from the northwest. 

The fifth week of the TOWARD experiment was from 29 
October to 5 November. Windspeeds were moderate on the 29th 
and 30th and were from the northwest. A wind shift and 
increase in windspeed began at approximately 0900 PDT on the 
3ist. Evaluation of Symoptic@™enharts fommeemoctober indicated 
the presence of a large-scale pressure gradient associated 
with an oncoming front. On November 2nd, the winds were 
again out of the northwest and San Diego returned to the 
typical sea breeze pattern. 

The development of a strong high pressure center 
over the western United States produced the low windspeeds 
and shift in direction evident on 4 and 5 £4November. 
Windspeeds ranged from approximately 2 to 6 m/s. 

Data collection for the final week of Phase I 
concluded on 7 November. Winds on the 5th were light and 
southwesterly. By the 7th, the typical sea breeze pattern 
again prevailed. Windspeeds were low and out of the west in 
early morning, but shifted to the northwest by early after- 


noon. A maximum value of about 5 m/s was recorded at 1105 
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Figure 3.6 
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Figure od Surface analyses 1800 GMT 23 (a) 
nd 25 (b)) October 17a 
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Figure 3.11 Surface analyses 1800 GMT 3 (a) 
and 4 (b) November 1984. 
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PDT followed by steadily declining values over the remainder 


of the sampling period. 
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2. Phase II 


Data collection for Phase II of the TOWARD experi- 
ment began on 12 March 1985. It can be seen from Fig. 3.14 
that winds on the 12th and 13th were light to moderate and 
were primarily westerly. On March 15th and 16th winds were 
again moderate, but had a more typical northwesterly flow. 

Synoptic charts for this period do not indicate any 
Significant weather activity. Conditions were dominated by 
an offshore high pressure center throughout the week. 

On Marea elcthes synoptic eharts depict a tight pres- 
Sure gradient resulting from the passage of a low pressure 
center through the San Diego area. The result was a sharp 
increase in windspeed. A maximum value of roughly 10 m/s was 
recorded at 1550 on that date. During the remainder of the 
week, the development of a strong offshore ridge reduced 
windspeeds and produced a northwesterly flow. 

The final week of the TOWARD experiment covered 25 
homo Mateus etacual wind shift accompanied by light to 
moderate winds can be seen on the 25th. Synoptic charts for 
that date show a front passing through the San Diego area 
between 1000 and 1600 local time, at an approximate rate of 
4-5 m/s. As it moved westward, the low deepened and the 
frontal band remained southeast of San Diego through 1000 on 
the 26th. The winds reflect this frontal movement. A shift 
in direction to a more typical northwesterly flow occurs in 
late morning and a sea breeze pattern is apparent by late 
afternoon. The maximum wind speed recorded was about 7 m/s 
age U7 LD. 

Winds on the 27th were out of the west-southwest and 
had generally higher magnitudes that those recorded the 
previous day. These values ranged from a mimimum of about 3 
m/s to a maximum of just over 8 m/s at 1015. This maximum is 
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large-scale pressure gradient associated with a deep low 
center over Vancouver, Canada. This system quickly moved 
southeast, where the low center deepened over Arizona and 
New Mexico for the next several days. The result was a tight 
pressure gradient which produced the high wind values 
collected on 28 March. On that date, a maximum wind value of 


over 12 m/s was recorded at 1415. 
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IV. DRAG COEFFICIENT RESULTS 


A principal motivation of this thesis is to examine the 
variation of the sea-surface drag coefficient with varying 
wind and wave states. Since windspeed magnitudes were gener- 
ally low to moderate (i.e. < 10-11 m/s) during both phases 
of TOWARD, the primary emphasis will be placed on Cpy varia- 
tions as a function of wave state. 

The question of whether the neutral drag coefficient is 
constant or increases with windspeed, wave state and fetch 
has received considerable Study in recent years. Garratt 
(1977) summarized sea-surface drag coefficients from 1/7 
experiments (based on eddy correlation and profile methods) 
published between 1967 andgi)/5. Obsermvatvons cf SC Greer. 
sentative of Vavhereht sof 10s ee loube daa tone: tonmot Ure. 
showed that Charnock's relation (Eqn. 26) was valid over a 
windspeed range of 4 to 21 m/s. 

Large and Pond (1981) compared measurements of the 
momentum flux obtained by both Reynolds flux and dissipation 
methods. They reported that a comparison of both methods 
Showed excellent agreement, on average, for windspeeds of 
4-20 m/s. Large and Pond also plotted neutral drag coeffi- 
clents against U;g. Although they found a slight increase in 
the average €7,, with windspecds@be lows Ogin Sycieys felts tian 
this relationship could be adequately described by a 


constant, where: 
10~G@ae : jie 4<Uj9<10 m/s (59) 


Above 10 m/s, the values of Cyn demonstrated a more rapid, 
linear increase up to windspeeds of 26 m/s. 
Smith (1980) used the eddy correlation method to obtain 


Windspeed and heat flux measurements over the windspeed 
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range of 6-22 m/s. He noted that the Charnock relation also 
held for neutral sea-surface drag coefficients measured for 
long fetch. These Cpy values were clearly seen to increase 
with increasing windspeed. Several points were not 
adequately explained by the Charnock hypothesis, however. 
These included: 

l. The dependence of the drag coefficient on windspeed 
was greater than predicted if the data at low windspeeds 
were not included. 

2. If the data in Garratt's summary are viewed sepa- 
rately at windspeeds both above and below 10 m/s, the drag 
coefficient also appears to be nearly constant at windspeeds 
below 10 m/s. Smith's average Cpy value for low windspeeds 
was (75a alee At higher windspeeds, this coeffcient tends 
to increase rapidly as a function of windspeed. 

Thus, all three studies indicate a relatively constant 
Cpn value at windspeeds below 10 m/s. It follows then that 
variations in the neutral drag coefficient at low to 
moderate windspeeds must be due to changes in the sea state. 

Garratt did not explore the behavior of Cpy with respect 
to fetch or sea state. He felt that current quantitative 
information was insufficient to adequately determine the 
dependence of the drag coefficient on factors other than 
windspeed. 

Large and Pond found the Cpy to be largely independent 
of fetch. For both long and short fetches, the mean drag 
coefficient remained at 1.14 x 10° for windspeeds less than 
HO m/s. 

Smith's comparison of Cpy with fetch showed higher drag 
coefficients at short fetches than those at longer fetches. 
This is to be expected since at short fetches actively 
growing waves absorb momentum from the wind while waves at 
longer fetches tend to be more in equilibrium with the wind 
Buea, 
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Smith also considered the variation of the neutral drag 
coefficient with™sea state. A linear eter ressiongOume ,  jeeae 
rms-wave height showed an increase in Cpy with an increase 
in wave height, although the correlation coefficient was 
rather low (.58). 

More recent studies have included a drag coefficient 
model discussed by Geernaert (1985a). Geernaert found the 
model exhibited considerable sensitivity to fetch. The 
influence of fetch on the model calculations was due to the 
wave spectral variations incurred by the dependence of the 
Phillips coefficient on the wave saturation parameter, 
Co/Us, where Co represents the phaSe speed. Geernaert 
observed that high drag coefficients often resulted from 
fetch limited condi tionse 

Wu (1985) also looked at fetch as influenced By the 
Phillips spectral coefficient. He found larger spectral 
coefficients at shorter fetches than those at longer 
fetches. 

Measurements of the neutral drag coefficient over the 
sea by all methods often yield considerable scatter and the 
data must be carefully edited. Several types of data have 
been removed from this analysis based on dissipation 
methods. They include: 

Ll. Data collected during rain or fog conditions. Water 
droplets striking the hot film result in latent heat 
release. This sudden cooling induces a large jump in 
voltage, producing spikes in the Wdataweetin. the inert 
Subrange particularly, noise from droplet contact contami- 
nates the high frequencies and results in very large drag 
coefficient estimates. 

2. Periods of light “winds also produce erratic” Cpa 
values as a result of convection off the tower mast. 

3. Flow distortion is a problem with light or moderate 
Winds, particularly from the Wsoubh seinem beste di receilome oe 


data collection is from west  touneweae 


62 


4. Radio transmissions, particularly in the autumn, 
also posed problems. The hot films serve as antennae for 
transmitted energy, occasionally producing unreliable Con 
values. 

Analysis was based on plotting Con as 2a function of Ujg, 
as in the earlier studies by Garratt, Large and Pond, and 
Smith (Fig. 4.1). The data were then separated into two sets 
in order to compare short and long fetch conditions. Short 
fetch waves occurred when winds originated from north to 
east (Fig. 4.2), while long fetch conditions were approxi- 
mated during west to northwesterly winds (Fig. 4.3). 


Ug 


TOWARD 





Figure 4.1 Con VS Ujg for long and short fetches. 


Con Was again plotted against Ujg for each fetch condi- 
tion. The mean and standard deviation for each 1 m/s 
velocity band was calculated and can be found in Table III. 

Table III indicates that higher drag coefficient values 
did result from short fetch conditions. Mean Cpy values for 
both the 2-3 and 3-4 m/s bands under short fetch conditions 


significantly exceeded those for long fetch. 
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TOWARD 


WIND DIRECTION 380-115 





Figure 4.2 Con VS Uyo for short fetch conditions. 


TOWARD 


‘NINO DIRECTION 250-330 





Figure 4.3 Con VS Uyg for long fetch conditions. 


Con Values for long fetch conditions also varied and 


were not constant as earlier predicted. Magnitudes were 
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TABLE IV 


Cyn vS. Ujg for Long and Short pet ene: 


SHORT FETCH CONDITIONS: 360-115 DEGREES 


me 0>2.0 m/s 2.0-3.0 m/s 3.0-4.0 m/s 
(1 point) (8 points) (2. “pemmts)) 
acs: Sy Aas 


S.d.=--- 


mean=2.017 
Shaw = .Sol 


mean=1.656 
Sea.-.407 


LONG FETCH CONDITIONS: 250-330 DEGREES 


1.0-2.0 m/s 

cl porn ) 

Benue- Jie NALS, 
Bic 


4.0-5.0 m/s 
(79 pointes ) 


2.0-3.0 m/s 
(se potmts ) 
mean=1.485 
s.d.=.828 


5.0-6.0 m/s 
(Geepolnes ) 


3.0-4.0 m/s 


Col pointes ) 
mean=1.432 
Sud. = G50 


bi = 770) mys 


(Zo. DOinte sy) 


mean=.955 ne =.980 mean=.896 
Sed amo Sia = 94 oe s.d.=.377 
7.0-8.0 m/s 8.0-9.0 m/s 9.0-10.0 m/s 
(14 points) Com OGHnit,S ) (4 points) 
mean=.856 mean=.939 mean=1.413 
s.d.=.463 s.d.=.254 Sac. —. 104 


10.0-11.0 m/s 
(ib eexe st ae) 


value=1.479 
S Grae 


fairly high for the 1-4 m/s bands, dropped off between 4-8 


m/S, and subsequently increased between 8-11 m/s. 
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Overall, the dependence of the drag coefficient on wind- 
Speed was found to be insignificant above 4 m/s, based on 
Scatter in the data as well as the large standard deviation 
at each windspeed interval. Below 4 m/s, each there was an 
apparent increase in the Cpy with a decrease in windspeed. 
These values are somewhat suspect, however, due to the 
uncertainty in the quality of the data at very low 
velocities. 

Daily time series of wind stress were also examined, 
with both Cpy and Uy plotted against time. Lack of data, low 
winds or strong stability prevented a close examination of 
the wind stress in some instances. Several case studies are 
described below, however. Appendix A contains figures rele- 


Vant Eor tne aisecussionm: 
1. Frontal Cases: 15, 16-18 October and’ 26-29 Maren 


Fig. A.6 shows the results when a front was passing 
the tower on 15 October7 9G, values (oanged ait mone. mec 
2.058. Magnitudes can be seen to rise with rising winds 
until 0900, fall with decreasing winds between 0900 to 1030, 
and remain fairly steady with the winds near constant 
values. A sharp drop in windspeed at 1200 did not produce a 
COrrespondime drop ip Con values. Cpn increased during this 
period of low winds to a maximum at 1500, and then dropped 
off sharply. Synoptically. Fie- A.G00shews (a) winds sha 
between 1400-1500 for southerly to westerly winds with a 
corresponding increase in windspeed. Con values began rising 
shortly afterward and continued to increase up to the end of 
the sampling period depite a drop in windspeed beginning at 
approximately 1700. 

A similar frontal case is presented in Fig. A./ for 
i628 Octoser. Both Us and Cyn can be clearly seen to 
increase with rising winds to a maximum at 0600 on the 18th. 


A wind shift from southerly to northwesterly winds occurred 
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Shortly after 6 am and several sharp drops in windspeed took 
- place before the wind leveled off at approximately 6 m/s. 
The drag coefficient shows a similar trend but does not 
reflect the rapid decreases evident in the windspeed between 
0500 to Q800. 

On 26-29 March, passage of a front as well as south- 
erly winds produce high Con values with considerable scatter 
ere weed . 3uly).. 

In each of these frontal cases, Cpn increased just 
before passage of the front. This trend is consistent with 
an earlier study by Denman and Miyake (1973) in which drag 
coefficients were seen to increase on the leading edge of a 
storm, and then decrease slightly or remain constant. They 
attributed this behavior to the nature of the wave field. 
Geernaert (1985b) also observed higher drag coefficients 
preceding an approaching front. He postulated that wave 
fields far ahead of and behind cold fronts are _ semi- 
independent of one another and are generally in equilibrim 
With local winds. Waves generated behind the front, often in 
a direction perpendicular to the baroclinic zone, can under 
certain circumstances travel through the front and interact 
with the wave field just ahead of it. Since the waves in the 
warm sector often run parallel to the local winds, the 
convergence of the two fields produces a much more energetic 
Sea state just ahead of the front. The result would be 


higher Cpyj estimates. 


femcea Breeze Cases: 18, 27, 29 October and 13, 15, 18 
March 


The sea breeze effect is apparent in much of the 
data from the TOWARD experiment. On each of the above dates 
Mites, Ani2yeAgl4, A.225e4.23 and Ag25),ethe Cpy rises 
Sharply and then decreases slowly to near equilibrium. 


Described in terms of surface roughness, drag coefficient 
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values are known to increase with growing wave state. As the 
wave field becomes increasingly saturated with momentum, 
absorption of energy slows, resulting in drag coefficients 


that approach a smaller and more constant value. 


3. Shift from Long to Short Fetch Conditions: 3-5 


November 


Data from 3-5 November depict a shift from long to 
short fetch conditions.) Long” feten seondtetons ee epee ome 
until approximately 2300 on 3 November, as indicated by 
northwesterly winds. The drag coefficient estimates (Fig. 
A.18) are correspondingly low. With the shift to northerly 
winds at 2300, a sharp increase can be seen in Cpy. Although 
these values decline over the next three hours, they again 
increase due to the variable wind direction produced by a 


Series Gf wind) Shiites  oMmmeneeaitn. 
4 Variable Winds a, hiewemp ens 


High Cpy values as result of “variable wand direc 
tions can be seen between 0900 to 1200 on 7 November (Fig. 
A.20). Magnitudes become more constant as winds shift to a 


Steadier direction. 


5. Long Fetch and a Con Decrease wlth Times lo.19, 22 
March 


These examples indicate long fetch conditions due to 
the northwesterly winds, as well as a general decrease in 
Con With time (Figs. A.24, A.26 and A.29). Since long fetch 
waves are known to be more steady state, they tend to 
absorb less momentum than would short fetch waves. 
Consequently, the drag coefficient estimates have lower 
magnitudes and tend to approach a constant value more 


rapidly than those seen in earlier examples. 
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In summary, case studies suggest a Cpy dependence on 
Sea state as illustrated by both surface roughness and the 
@eecee OL equilibrium between Uj, and the waves (i.e. the 
amount of momentum entering the wave field for growth). 
Basic to this hypothesis on wind-wave coupling is that Cpy 
is generally independent of Ujg for windspeeds between 4-9 
m/s, and that the large amount of scatter is associated with 
systematic trends in sea state, such as frontal and sea 
breeze conditions. This scatter should be reduced once wave 
State is incorporated as an additional parameter in Cpyj 


Goanme lations . 
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V. SUMMARY AND CONCLUSIONS 


Surface layer measurements covering a wide range of 
windpeeds have been analyzed for both phases of the TOWARD 
experiment. These measurements have been correlated with 
relevant mesoscale and synoptic features. The results were 
then interpreted relativestosthe Gn. 

Neutral drag coefficient values for both phases of 
TOWARD showed no significant windspeed dependence in the 
range between 4 and 9 m/s. The mean value was determined to 
be (.94 + 0.4)10°, which is somewhat smaller than results 
obtained by Smith (1980) and Large and Pond (1981). These 
differences appear to be insignificant, however, due to the 
scatter present in all three “data Ssebs Seay ote s ee 
higher for windspeeds less than 4 m/sec, but these values 
are suspect as low windspeeds often produce unreliable drag 
coefficient estimates. 

Definite trends were observed in Gay with fet cheanaeeea 
state. Short fetch conditions produced generally higher Cpy 
values than those for long fetch. Trends in sea state such 
as frontal and sea breeze conditions were reflected in drag 
coefficient estimates as well. As with Large and Pond, Cpy 
estimates were seen to increase sharply just prior to the 
passage of a front. During sea breeze Neonditlous rice. 
values increased with time. As windspeeds and the wave field 
reached equilibrium, however, the drag coefficients 


decreased toward amore constant value. 
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APPENDIX A 
TOWARD DAILY TIME SERIES 


This appendix contains daily time series for wind speed 
and direction, air and sea surface temperatures, the neutral 
drag coefficient and Us, for both phases of the TOWARD 
experiment. These figures provide a reference for the 


discussion provided in chapters three and four. 
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